A single defect approximation for localized states on random lattices 
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Geometrical disorder is present in many physical situations giving rise to eigenvalue problems. 
0^ ' The simplest case of diffusion on a random lattice with fluctuating site connectivities is studied 

0^ I analytically and by exact numerical diagonalizations. Localization of eigenmodes is shown to be 

CTS . induced by geometrical defects, that is sites with abnormally low or large connectivities. We expose 

a "single defect approximation" (SDA) scheme founded on this mechanism that provides an accurate 
quantitative description of both extended and localized regions of the spectrum. We then present 
a systematic diagrammatic expansion allowing to use SDA for finite-dimensional problems, e.g. to 
determine the localized harmonic modes of amorphous media. 
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fH I Since Anderson's fundamental work W , physical systems in presence of disorder are well known to exhibit localiza- 

Q tion effects Q. While most attention has been paid so far to Hamiltonians with random potentials (e.g. stemming 
1^ from impurities), there are situations in which disorder also originates from geometry. 

H , Among these are of particular interest the harmonic vibrations of amorphous materials as liquids, colloids, glasses, 

' ' ... around random particle configurations. Recent experiments on sound propagation in granular media fQJI^ have 
C^ i stressed the possible presence of localization effects, highly correlated with the microscopic structure of the sample. 
"^ ' The existing theoretical framework for calculating the density of harmonic modes in amorphous systems was developed 
I t ' in liquid theory. In this context, microscopic configurations are not frozen but instantaneous normal modes (INM) 
Cd ' give access to short time dynamics [|| . Wu and Loring g and Wan and Stratt 0] have calculated good estimates of the 
H ' density of INM for Lennard- Jones liquids, averaged over instantaneous particle configurations. However, localization- 
I delocalization properties of the eigenvectors have not been considered. 

Diffusion on random lattices is another problem where geometrical randomness plays a crucial role |q]. Long time 
dynamics is deeply related to the small eigenvalues of the Laplacian on the lattice and therefore to its spectral 
fj ' dimension. Campbell suggested that diffusion on a random lattice could also mimic the dynamics taking place in a 
complicated phase space, e.g. for glassy systems M. From this point of view, sites on the lattice represent microscopic 
configurations and edges allowed moves from a configuration to another. At low temperatures, most edges correspond 
to very improbable jumps and may be erased. The tail of density of states of Laplacian on random graphs was studied 
j^ , by means of heuristic arguments by Bray and Rodgers [|lO[. Localized eigenvectors, closely related to metastable 
^f^ ' states are of particular relevance for asymptotic dynamics. 

f^ , Remarkably, the above examples lead to the study of the spectral properties of random symmetric matrices W 

CN ' sharing common features. In amorphous media, the elastic energy is a quadratic function of the displacements of the 
^^ ] particles from their instantaneous "frozen" positions. The INM are the eigenmodes of the stiffness matrix W. As for 
^Zf^ ' diffusion on random lattices, W simply equals the Laplacian operator. In both cases, each row of W is comprised 
— .^ of a small (with respect to the size N of the matrix) and random number of non-zero coefficients Wij and most 
"jrt importantly, diagonal elements fiuctuate : Wu = — J2j(M) ^ij 10 • 

ri _ In this letter, we present a quantitative approach to explain the spectral properties and especially localization 

I . effects of such a random matrix W in the simplest case, that is when all off-diagonal elements of W are independent 

''O ' random variables. Our analytical approximation is corroborated by exact numerical diagonalizations. We then expose 

Ch , a systematic diagrammatic expansion allowing for the study of more realistic models in presence of correlated Wij 's. 

O ■ The spectral properties of W can be obtained through the knowlegde of the resolvent G{\ + ie), that is the trace 

of ((A + ie)l — W)~^ Ig]. Denoting the average over disorder by (•), the mean density of states reads 
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p(A) = lim Im G'(A + ie) . (1) 

The averaged resolvent is then written as the propagator of a replicated Gaussian field theory [g| 
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where Zi = z{4>i) = exp ( ^(A + ie)(j)^ 
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Replicated fields <j)i are n-dimensional vector fields attached to each site i. To ligthen notations, we have restricted in 

(^ to the scalar case. We shall focus later on W having an internal dimension, as happens to be in the INM problem. 

In the uncorrelated case, the Wij's (i < j) are independently drawn from a probability law V . To take into account 



geometrical randomness only, we focus on distribution V (Wij) = (l 
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Such a 



bimodal law merely defines a random graph : i and j can be said to be connected if and only if Wij does not vanish. 
Due to the scaling of the edge- probability -^, the mean site-connectivity q remains finite for large sizes N. We rescale 
the eigenvalues by choosing w — —^ to ensure that the support of the spectrum is positive and bounded when q ^ oo 




FIG. 1. Density of states p{\), inverse participation ratio w'^{X) and connectivity of the centers c(A) (divided by g) averaged 
over 2000 samples for q = 20, A'^ = 800 (all solid lines). The bottom parts of the oscillations of w'^ at small and large A 
are suppressed due to the statistical shortfall of eigenvalues. Dotted and dashed lines respectively show the EM A and (the 
extended region of) the SDA spectra. Inset : oscillations of w*(A) for A^ = 100, 200, 400, 800, 1600 and 3200 from top to down 
and fluctuations of c'(A) (divided by q) for TV = 100, 200, 1600 (same symbols) and 0.4 < A < 0.5. 



Numerical diagonalizations of the random Laplacian W have been carried out for different sizes, up to iV = 3200. To 
each eigenvector tpi^g of eigenvalue A^ normalized to unity is associated the inverse participation ratio wf = ^^ |V'i,f 1^- 
We then define u;*(A)dA as the sum of wj over all ipi^i lying in the range A < A^ < A -t- dX, divided by the number 
Np{X)dX of such eigenvectors. Fig. 1 displays p{X) and w^(A) for a mean connectivity q — 20 |13). The central part of 
the spectrum (A_ < A < A+) has a smooth bell shape and corresponds to extended states. For increasing sizes N and 
at fixed A, w^ vanishes as 1/iV and the breakdown of this scaling identifies the mobility edges : A_ ~ 0.47 ± 0.01 and 
A+ ~ 1.67 ± 0.03. Outside the central region, that is for small or large eigenvalues, the eigenstates become localized 
and the density exhibits successive regular peaks. 

We have measured for each eigenvector ^i^i the connectivity ci of its center, that is the site iq with maximum 
component |V'io,^l- The mean connectivity c(A) of the centers of eigenvectors having eigenvalue A is plotted fig. 1. It 
is a smooth monotonous function of A in the extended part of the spectrum. In the localized region, c(A) is constant 
over a given peak and integer- valued {c < C- = 10 on the left side of the spectrum, c > c+ = 33 on the right side); 
the center connectivity abruptly jumps when A crosses the borders between peaks. Furthermore, Table 1 shows the 



good agreement between the weight of peak associated to connectivity c and the fraction of sites having c neighbors, 
given by a Poisson law of parameter q |14| . 
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TABLE I. Weights p {i.e. integrated density of eigenvalues belonging) of the peaks and corresponding eigenvalues A obtained 
from numerical simulations (NUM) and SDA for low connectivities c. \^'^^' is measured at the top of the peak with absolute 
error ±0.0025 whereas the relative error on p^^^' is about 10% (except for the c = 10 peak). A similar agreement is reached 
for large connectivities. 



Therefore, numerics indicate that localized eigenvectors are centered on geometrical defects, that is on sites whose 
number of neighbors is much smaller or much larger than the average connectivity. To support this observation, it is 
intructive to consider a simpler model including a unique defect i. e. a Cayley tree with connectivity c for the central 
site and q + 1 for all other points |lj] (locally a random graph is equivalent to a tree since no loops of finite length 
are present). Looking for a localized state ipt with a radial symmetry -tpi = ipd{i) where d{i) is the distance between 
site i and center, the eigenvalue equations read c{ipo — "01) = g A V'o and (q + l)ipd ~ q i^d+i ~ V-'d-i = q ^ 'ipd for 
d> 1 [|5| . The eigenvalue problem reduces to the search of the solution of a homogeneous linear difference equation 
of order two (eqn. for d > 1) fulfilling a boundary condition (eqn. for d = 0). After a little algebra we have found 
that strong defects, such that |g — c| > y/q give rise to localized states around the central site with an eigenvalue 

A = -(1 ^). The predicted connectivities at mobility edges (c_ = 15 and c+ = 25 for q = 20) are in poor 

agreement with numerical findings. A more refined picture requires to take into account the connectivity fluctuations 
of the neighbours of the central site. We have thus considered a Cayley tree with a coordination number c for the 
central site, c' for the nearest neighbours and q + 1 for all other points. We have found that localized states due to weak 
(respectively strong) central connectivity c can disappear and become extended if the connectivity of the neighbors 
c' reaches large (resp. small) values. In other words, a defect can be screened by an opposite connectivity fluctuation 
of its surrounding neighbors. Numerics supports this scenario. As A varies, w'^{X) exhibit oscillations (interpreted 
as finite-size contributions coming from extended states) of rapidly decreasing amplitudes with increasing N. These 
oscillations are correlated (positively for small c and negatively for large c) with the fluctuations of the neighbors 
connectivity c'(A) around its mean value g + 1, see inset of flg. 1. 

Let us see how the above results may be recovered from theory. Due to the statistical independence of the Wij 's, the 
Uij (^ interactions are averaged out separately |10| . The resulting theory is of course invariant under any relabelling 

of the sites i and depends on the flelds (pi through the density p((/)) of sites i carrying flelds (fi = (f only |16|. The 
functional order parameter /9(0) is found when optimizing the "free-energy" |10| , [16[ 
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under the normalization constraint / d(pp{4>) = 1; z{(j)) has been defined in (0). This order parameter is simply related 
to the original random matrix problem through 
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The Ci are normalization constants going to unity as n vanishes. Therefore, the averaged resolvent reads G{\ + ie) = 

-i lim„^o / rf^p(^) (0^)^- 

Finding an exact solution to the maximization equation 5 In 'E./Sp{4>) = seems to be a hopeless task. This is a 
general situation, which arises in the study of the physics of dilute systems (for the case of sparse random matrices 
see, for example, 0,0 )• 



Identity (||) may however be used as a starting point for an effective medium approximation (EMA) . In the extended 
part of the spectrum, we expect all matrix elements appearing in (Q) to be of the same order of magnitude and thus 
p{(f>) to be roughly Gaussian. EMA is therefore implemented by inserting the Gaussian Ansatz 

p^M\$) = [2^^g{X))-'^e^^(^^\ , (6) 

into functional S (Q). The average EMA resolvent g is then obtained through optimization of lnS[(7(A)]. The resulting 
spectrum, which is given by the imaginary part of g divided by minus tt, is shown fig. 1. As expected, EMA gives 
a sensible estimate of the spectral properties in the extended region and of the mobility edges A^*^'^ — 0.468, 
\EMA _ 2^ 732. However, EMA is intrinsically unable to reflect geometry fluctuations and thus the presence of 
localized states Q. 

To do so, we start by writing the cxtremization condition of InS over p as 



p{ci>) = n[pm , (7) 



where the functional TC may be expanded as 
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h is a. multiplicative factor equal to unity in the n — > limit. Equations MM) describe an elementary lattice of one 
central site connected to k neighbors according to a Poisson distribution of mean q. Neighbors carry information 
about the random matrix elements through the order parameter p{ip) m) and interact with the central site via kernel 
exp(— 1((/) — ip)^/2q) (13,0). Self-consistency requires that the resulting order parameter at central site equals p [18|. 
Baring in mind the localization mechanism unveiled in previous paragraphs, we propose a single defect approximation 
(SDA). SDA amounts to make the central site interact with k neighbours belonging to the effective medium defined 
above. Since EMA precisely washes out any local geometrical fluctuation, we partially reintroduce them by allowing 
the connectivity k of the central site (the defect) to vary. The SDA order parameter is thus obtained through an 
iteration of equation (|7|) 

p^^^i$) = n[p''''^m . (9) 

Using the EMA resolvent g (0), we have computed the SDA spectrum for q = 20. The SDA extended part is shown to 
be in better agreement with numerical results than EMA on fig. 1 . Improvement is even more spectacular for localized 
states that were absent within EMA. We have found Dirac peaks whose weights and eigenvalues are listed Table 1. 
The agreement with numerical results is quite good. We have verified analytically that SDA peaks do correspond to 
localized states by calculating lime^o[e G(A + ie)G{X — ie)] |1^] using SDA with two groups of replicas. This quantity 
gives also access to w'^{X), whose value seems slightly higher than numerical measures close to the mobility edges. 

Starting from any sensible p, successive iterations of equation ([7|) would also provide more and more accurate 
descriptions of the "fractal" structure of the localized peaks at a price of heavier and heavier calculations. Besides 
being theoretically founded, SDA has the advantage that a single iteration from EMA (which is easily computable) 
succeeds in capturing localized states in a quantitative way p^ . 

In general, the components of W are correlated and the average over disorder requires an expansion in terms of 
the connected correlation functions of the Uij interactions 



Ylil + uy) = exp ^ u^j + 2 XI "y "fe''' + ■ • • ' (10) 

where the prime indicates that the sum runs over different pairs of sites. Free-energy (0) corresponds to the case 
where all terms in (|l^) but the first one vanish. The presence of these cumulants (of order 2,3,... in u) will result 
in the addition of cubic, quartic, ... p{(j)) interactions terms to InS. Though calculations become more difficult, the 
existence of a variational free-energy In S is preserved. This is all what is needed to derive EMA and the optimization 
equation (||). 

Let us see how SDA can be implemented to determine the INM spectrum of amorphous media. We shall restrict 
to liquids but our approach could also be applied to glasses using the formalism recently developed by Mezard and 



Parisi |2^. Particles i are individuated by their positions x^ and interact through a two-body potential F(xi — Xj) 
(hereafter bold letters will denote vectors in the D-dimensional real space). For a given microscopic configuration, 
INM are the eigenmodes of the D x A^-dimensional matrix Wij = d^V{Ki — Xj) / dxidxj {i ^ j). The calculation 
of the spectrum and the average over particle configurations (with the equilibrium Boltzmann measure at inverse 
temperature /?) can be performed at the same time by introducing a generalized liquid [pp] . Each particle is assigned 
a "position" r^ = (xi,0i) and the generalized fugacity reads z*{ri) = y z{(j)i) where y is the liquid fugacity and z 
is defined in (0). The grand-canonical partition function S may then diagrammatically expanded in powers of the 
Mayer bond b{ri, Vj) — exp(— /3]/(xi — Xj) -|- ^Wij{(j)i — 4>j)'^) — 1 (the summation over the D"^ internal indices of W is 
not written explicitely for sake of simplicity). With these notations, S coincides with formula (3.21) of Ref. |£l|. It is 
now straigthforward to take advantage of the variational formulation of the diagrammatic virial expansion by Morita 
and Hiroike |£1[. The generalized density p(r) — p{ii., 4>) of particles optimizes 



lnS[p]= j dvp{v)[\nz*{r)-\np{v) + l]+S , (11) 



where S is the sum of all diagrams composed of bonds 6(r, r') and vertices weighted with p(r) that cannot be split 
under the removal of a single vertex, see equation (4.6) of Ref. [pi|. Due to translational invariance in real space, p 
does not depends on x and we are left with a variational functional S of the density p{4>). Note that ( [ll|) contains (Q) 
as a special case when S includes only the simplest single bond diagram. The random graph model we have studied 
in this letter may be seen as a physical system for which keeping the first coefficient of the virial expansion only is 
exact. 

To our knowledge, Morita and Hiroike's work has not been used so far in the context of INM theory as a short cut to 
avoid tedious diagrammatical calculations. In addition, the variational formulation of |2l| ] allows to implement SDA 
in a practical way. We are currently attempting to apply the present formalism to characterize localized eigenstates 
in two- and three-dimensional granular media. 
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